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ELECTROMAGNETIC  TRANSMISSION  THROUGH  A 


SLOT  IN  A  PERFECTLY  CONDUCTING  PLANE 

I.  INTRODUCTION 

Although  the  computer  program  iu  [1]  can  treat  the  case  where  the 
rectangular  aperture  is  one  subsection  wide,  i.e.,  L^  •  1  in  Fig.  1  of  [1] 
this  computer  program  is  long  and  complicated.  A  computer  program  was 
written  specifically  for  the  L^  ■  1  case.  Relatively  short,  this  program 
is  described  and  listed  here.  It  consists  of  a  main  program  and  the  sub¬ 
routines  YMAT,  PLANE,  DECOMP,  and  SOLVE. 

The  formulas  that  are  programmed  are  presented  in  Section  II.  The 
main  program  is  described  and  listed  in  Section  III,  the  subroutine  YMAT 
in  Section  IV,  the  subroutine  PLANE  in  Section  V,  and  the  subroutines 
DECOMP  and  SOLVE  in  Section  VI. 


II.  FORMULATION 


The  magnetic  current  M  on  the  z  <  0  side  of  the  aperture  (see  [2, 
Fig.  2])  is  expressed  as 


L  -1 
x 


M  -  I  V  £ 

j-1  2  2 


(1) 


where  {V^}  are  unknown  coefficients  and  {M^}  are  expansion  functions 
defined  by  [1,  Eq.  (10)] 

M3;  -  u  T.(x)  P(y)  ,  J-1.2....L-1 

— j  “*  J  * 


(2) 


Here,  u  is  the  unit  vector  in  the  x-directlon,  T.(x)  is  the  triangle 
-x  j 

function  defined  by 


When  the  Incident  magnetic  field  is  given  by  (5)»  the  coefficients 


{Vj}  in  (1)  are  the  elements  of  the  column  vector  V  that  satisfies 


YV  -  -  (PlnC) 


9y 


(7) 


where  the  ith  element  of  the  column  vector  (Pinc)py  is  called  (p^QC)gy 
and  is  given  by  [1,  Eq.  (59)] 

.  2 

.  .  -inc 

.  kAx  cos  9  . 

(Pjnc)fi  -  2AxAy  sin  9inc  (~ - f - )  .***  cos  9  . 

1  07  kAx  cos  9^ 

2 

i-l,2,...Lx-l  (8) 

In  (7),  Y  is  the  square  matrix  whose  ij  element  is  called  and  is  given 
by  [1,  Eq.  (23)] 

Yij  "  [I  "  I  Ix<J-i+1>  +  ^~1-23/2^  IcU-i+D 

+  \  I( j-i-1)  -  (^^-3/-2>  I  (j-i-1)  +  1  -x-  (I  (j-i+1) 

1  X  (kAx;  c 


-  2Ic(j-i)  +  Ic (j-i-1))] 


(9' 


The  j  that  appears  in  the  factor  jAxAy/Crrn)  on  the  right-hand  side  of 

(9)  is  J-l.  Each  of  the  rest  of  the  j's  on  the  right-hand  side  of  (9) 

is  the  subscript  j  on  Y^ .  In  (9),  n  is  the  impedance  of  free  space. 

Moreover,  I  and  I  are  given  by  [1,  Eqs.  (27)  and  (28)] 
c  x 


/x2+y2 


X 

u 

f 


(i  +  l/2)kAx 


x 

u 

x£  -  (i  -  l/2)kAx 

It  is  evident  that  I  (i)  is  even  in  i  and  that  I  (i)  is  odd  in  i. 

c  x 

When  the  incident  magnetic  field  is  given  by  (6) ,  the  coef¬ 
ficients  {Vj}  in  (1)  are  the  elements  of  the  column  vector  V  that 
satisfies 


YV  -  -  (?lnC)xx 

where  the  ith  element  of  the  column  vector 
and  is  given  by  [1,  Eq.  (65)] 


(P1"0)  is  called  (pj“c) 
'xx  i 


(pfc> 

i  xx 


-  2AxAy  (• 


sin  «»  ♦ 


inc 


kAy  cos  <t> 


Inc 


J  (kAy/2)  cos  d)1*10 
i-1.2,...Lx-l 


In  (9), 


2-Lx  <  j-i  <  Lx-2 
but,  because  (9)  is  even  in  j-i, 

0  <  j-i  <  Lx-2 


is  sufficient.  Therefore, 


is  sufficient  in  (10)  and  (11). 


The  integrals  (10)  and  (11)  are  evaluated  by  using  the  following 
four  term  approximation  [1,  Eq.  (39)] 


lT  =  e  ri[l-j(r-r1)  -  f  (r-r±)2  +|  (r-^)3] 


where 


l2~2 

r  -  Jx  +  y 


rt-  | ikAx 


Substitution  of  (20)  into  (10)  gives  {1,  Eq.  (42)] 


y«  -  <cini  +  C2D2  +  C3°3  +  W2e 


where 


u!  "  1  -  T  +  jri  (l  ■  T> 


u2  *  ri  -  j 


U3  -  -  f  (1  +  ir±) 

U4-6 


y  x 
u  u 

dy 

0  x„ 


dxr  t  n  ■  1, 2,3,4 


Substitution  of  (20)  into  (11)  gives 


i  (i)  -  (XjO  +  x2o2  +  x3u3  +  w  -T3T 


where  U^,  U^,  U3>  and  are  given  by  (24) -(27)  and 


X 

n 


dx  xr 


n-2 


n=l ,2,3,4 


Using  the  indefinite  integrals  [1,  Eqs.  (46)-(54)],  we  obtain 
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T  (xur4 
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<V4  '  Vj> 


*1  "  1  ‘VVV 
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u  xu  Z  xj. 
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r^y  kAx 
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,  (kAx)2  V 
4  3  ' 


sign  (i) 


(39) 

(40) 

- 

'*■ 

f41')  m 


where  sign  (i)  denotes  the  algebraic  sign  of  i.  If  i*0,  then  sign  (i) 
is  inconsequential  because  the  term  that  multiplies  it  is  zero.  In 
(31)-(38) > 


<7T 


■*) 


x  +  r. 


v  log  (j a - i) 

“  S  +  r3 


0 


(< 


Here,  log  denotes  the  natural  logarithm. 

Substituting  Y/2  for  in  [2,  Eq.  (27)],  we  obtain  for  the  com¬ 
plex  power  P  transmitted  through  the  aperture 


Pt-|v[YV]*  (■ 

where  V  is  the  transpose  of  V  and  *  denotes  the  complex  conjugate. 
Since  the  incident  magnetic  field  is  given  by  either  (5)  or  (6),  YV 
is  given  by  either  (7)  or  (15)  so  that  (44)  reduces  to 

Pt  “  “  I  VP*  0 


-►  +inc 

where  P  is  (P  )gy  if  the  incident  magnetic  field  is  given  by  (5)  and 
P  is  (P^ne)  if  the  incident  magnetic  field  is  given  by  (6). 

The  transmission  coefficient  T  is  the  ratio  of  the  real  power 
transmitted  through  the  aperture  to  the  real  power  P^nc  incident  on 
the  aperture.  Since  the  incident  magnetic  field  is  given  by  either 
(5)  or  (6),  we  obtain 


P.  *  HL  AxAy  cos  a  0 

inc  x 

where  a  is  the  angle  between  the  propagation  vector  of  the  incident 

wave  and  u  .  Here,  u  is  the  unit  vector  in  the  z  direction.  Since 
— z  — z 

the  incident  magnetic  field  is  given  by  (5)  or  (6),  inspection  of 
[1,  Fig.  1]  reveals  that 


cos  a  ■  -  sin  8 


( 


where  6  is  either  the  angle  6 
(45)-(47),  we  obtain 


in  (5)  or  the  angle  in  (6).  From 


Real (VP*) 
2nLxAxAy  sin  6 


(4J 


The  scattering  cross  section  r(6,<t>)  is  the  area  that  the  incident 
power  per  unit  area  must  be  multiplied  by  in  order  to  obtain  the  power 
which,  when  radiated  omnidirectionally  into  the  z  >  0  half  space,  would 
produce  the  actual  power  per  unit  area  at  (6,$).  The  above  definition 
of  x  leads  to  [1,  Eq.  (72)] 

,4 


t/X 


3  2 
327 rJn 


(4< 


*  "“ID 

where  X  is  the  wavelength  and  F  is  the  transpose  of  a  measurement 

vector  P®.  For  the  9  polarized  pattern  in  the  y»0  plane,  T  is  called 

(t)^  and  ?®  is  (P™)  ^  where  (P®)^  is  the  column  vector  whose  ith  ele- 

ixic 

ment  is  given  by  (8)  with  the  angle  of  incidence  9  replaced  by  the 
observation  or  "measurement"  angle  0®, 


<*VX 


k' 


3  2 
32tt  n 


(5( 


For  the  x  polarized  pattern  in  the  x«0  plane,  T  is  called  (T)xx  and 
P®  is  (P®)  where  (P®)  is  the  column  vector  whose  ith  element  is 

XX  XX 

inc 

given  by  (16)  with  the  angle  of  incidence  <{>  replaced  by  the 
measurement  angle  <£®. 


(x)  !\ 


k" 


(?\ 


Si  2 


7  9 


(53 


III.  THE  MAIN  PROGRAM 


The  main  program  uses  the  subroutines  YMAT,  PLANE,  DECOMP,  and 
SOLVE  to  calculate  the  coefficients  {v^ }  appearing  in  expression  (1) 
for  the  magnetic  current,  the  transmission  coefficient  (48),  and  the 
scattering  cross  sections  per  square  wavelength  (SO)  and  (51).  The 
main  program  is  described  and  listed  in  this  section.  Sample  Input 
and  output  data  are  provided  so  that  the  user  can  verify  that  the 
program  is  running  properly. 

In  the  listing  of  the  main  program,  line  4  defines  the  input 
data  file  and  line  S  defines  the  output  data  file.  The  input  data 
are  read  according  to  lines  10  and  11  which  are 

READ(20,11)  LX,  LI,  NTH,  DX,  DY,  TH 
11  FORMAT (313,  3E14.7) 

Here,  LX  is  L^,  DX  is  Ax/X,  and  DY  is  Ay/X  where  X  is  the  wavelength, 
and,  as  in  Section  II,  L^Ax  is  the  length  of  the  aperture  in  the  x 
direction  and  Ay  is  the  width  of  the  aperture  in  the  y  direction.  LX 
is  a  positive  integer  greater  than  or  equal  to  2.  LI  is  either  1  or  2. 
If  LI  is  1,  then  the  incident  magnetic  field  is  given  by  (5)  and  TH  is 
0  of  (5).  If  LI  is  2,  then  the  incident  magnetic  field  is  given  by 

inC 

(6)  and  TH  is  4)  of  (6).  The  input  variable  TH  is  in  degrees.  The 
normalized  cross  section  (50)  is  calculated  at 

6m  -  (J-l)ir/(NTH-l),  J  -  1,2,...  NTH  (52 

and  is  written  on  the  output  data  file  under  the  heading  TAUl.  The 
normalized  cross  section  (51)  is  calculated  at 


4>  -  (J-l)x/ (NTH-1) 


(53) 


9  J  "  1>»2|  •  •  •  9 NTH 

and  is  written  on  the  output  data  file  under  the  heading  TAU2.  The 
right-hand  sides  of  both  (52)  and  (53)  are  in  radians. 

Minimum  allocations  are  given  by 

COMPLEX  Y(N*N),  P(2*N),  B(N) ,  V(N) 

DIMENSION  IPS(N) 

where 

N  -  LX  -  1  (54) 

Line  20  stores  by  columns  in  Y  the  elements  of  the  matrix  yr^r  Y 

jAxAy 

where  Y  appears  in  (7).  Line  24  stores  in  P(l)  to  P(N)  the  elements  of 
_i_  (Pinc)Qy  where  (Pinc)Qy  appears  in  (7).  Here,  N  is  given  by  (54). 
Line  24  also  stores  in  P(N+1)  to  P(2*N)  the  elements  of  (P1®0)^ 

where  (pinc)^  appears  in  (15) .  Equation  (7)  is  recast  as 

'3Sk  Yl?  ’  J2,m  l2SE7  (?1”' )9yI  <55) 

Equation  (15)  is  recast  as 

YlY  •  ^2lTT1  ‘aEZF  <56> 

The  square  matrix  [  Y]  is  common  to  the  left-hand  sides  of  both 

(55)  and  (56) .  This  is  the  matrix  that  resides  in  the  computer  program 
variable  Y.  The  bracketed  quantity  on  the  right-hand  side  of  (55)  is 
the  column  vector  that  resides  in  the  computer  program  variables  P(l) 
to  P(N)  where  N  is  given  by  (54).  The  bracketed  quantity  on  the  right- 
hand  side  of  (56)  is  the  column  vector  that  resides  in  the  computer 
program  variables  P(N+1)  to  P(2*N).  If  LI  is  1,  DO  loop  16  performs 


the  multiplication  by  the  factor  UV  -  2jim  in  (55)  in  order  to  store 
the  elements  of  the  right-hand  side  of  (55)  in  B(l)  to  B(N).  If  LI 
is  2,  DO  loop  16  stores  the  elements  of  the  right-hand  side  of  (56) 
in  B(l)  to  B(N) .  If  LI  is  1,  lines  33  and  34  put  in  V(l)  to  V(N)  the 
elements  of  the  column  vector  V  that  satisfies  (55).  If  LI  is  2, 
lines  33  and  34  put  in  V(l)  to  V(N)  the  elements  of  the  column  vector 
V  that  satisfies  (56).  Since  the  elements  of  are  the  {V^}  of  (1),  it 
is  evident  that  the  coefficients  (V^ }  in  the  expansion  (1)  for  the  mag¬ 
netic  current  M  will  reside  in  V(l)  to  V(N) . 

Equation  (48)  is  recast  as 


Real  (V 


t  _ _ zAxfo 

ilL  sin  8 
x 


P*) 


(57) 


In  (57),  B  is  TH  and 


1  -+ 

P  is  the  column  vector  that  stored  in  either 


2AxAy 

P(l)  to  P(N)  or  P(N+1)  to  P(2*N)  according  as  LI  is  either  1  or  2.  With 


regard  to  (57),  DO  loop  17  accumulates  V 


2AxAy 

of  (57)  in  the  computer  program  variable  T. 


P*  in  U1 .  Line  44  puts  T 


Equation  (50)  is  recast  as 


<TVX 


(k2AxAy) 2 


3  2 
8"  TT 


— i —  (pm)  ^ 
2AxAy  '  '9y 


(58) 


Equation  (51)  is  recast  as 


(T)  /x2  .  OsJaM 


XX 


3  2 

8xJn 


— - —  (Pm)  v 

2AxAy  K  ;xx 


The  index  J  of  DO  loop  19  obtains 


(59) 


9®  -  (J-1)7t/(NTH-1) 


(60) 


in  (58)  and 


<j>®  -  (J-l)ir/ (NTH-1) 


(61) 


in  (59).  The  value  of  Che  right-hand  side  of  (61)  is  the  same  as  that 

of  (60).  Inside  DO  loop  19 ,  line  53  puts  this  common  value  in  TH.  With 

regard  to  (58)  and  (59),  line  54  puts  the  elements  of  yjgi&y  @y  in 

P(l)  to  P(N)  and  the  elements  of  (P®)^  in  P(N+1)  to  P(2*N).  DO 

loop  21  lies  inside  DO  loop  20  whose  index  is  K.  If  R  ■  1,  DO  loop  21 

accumulates  Y^Ay  1x1  111  •  If  K  ■  2,  DO  loop  21  accumilates 

—r-  (P®)  V  in  01.  When  K  ■  1,  line  64  puts  (t)ft  /X2  of  (58)  with  0m 
zcxuy  xx  ey  - 

given  by  (60)  in  TAD(l).  When  K  ■  2,  line  64  puts  (t)^/!2  of  (59)  with 


<j>  given  by  (61)  in  TAU(2). 
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00 1C  LISTING  OF  THE  flAIN  PRGG5AB 

002  COMPLEX  U,UV ,T  (1600) , P  (100) , B (40) #7(40) , 01 ,CCNJG 

003  DIHEISIOB  IPS  (40)  ,TAO  (2) 

004  CPEN  |UMIT*20,FILE* *  HAUT23. DAT*) 

00 S  OPED (OHIT-2 1 ,FILE- *  BAUTZ4. DAT* ) 

006  PI*3- 141593 

007  FT  A*  376*730 

003  0=  (0. ,  1.) 

009  07«2.*PI*BTA*0 

010  HEAD  (20# 1 1)  IX, LI, NTH ,  DX,DY,TH 

011  11  FORMAT (31 3, 3E1 4- 7) 

012  9BITE (21,12)  LX, II, NTH, DX, DY,TB 

013  12  FORIIAT  (*  LX  LI  HTH*  ,5X,  *  DX  •  ,  12X,  *  DY»,  12X,  *TH»/1X, 

014  1  313,3214.7) 

015  EK*2.*PI 

016  DX*DX*BK 

017  EY*DY*BK 

018  P8*180-/PI 

019  TH»TH/P8 

020  CALL  THAT  (LX,DX,DY,Y) 

021  HRITE(21,13)  (YU), 1*1, 3) 

022  13  FORMAT  (•  Y»/ (1X,€E11.4) ) 

023  N=LX- 1 

024  CALL  FLAKE (TH, LX,DX,DT ,P) 

025  KBIT  E  (2 1,14)  (P  (I),  I*  1,3) 

026  14  FORM XT {*  P»/1X,6I1U4) 

027  IA*(LI-1)*N 

02  8  IB*IA 

029  CO  16  J»1.H 

030  IB=IE+1 

031  E { J)  *UV*P  (IE) 

C32  16  CONTINUE 

C33  CALL  DECOMP  (N,  IPS,  T) 

034  CALL  SCL7E  (N,IPS,  Y,B,7) 

035  BRITE(21,24)  (V(I),I*1,N) 

C36  24  FOBH AT  I*  COEFFICIENTS  7  CF  N16NITIC  CORSINT  *, 

037  1  • EXPANSION  FDNCTIONS*/{ IX, 6111-4) ) 

C38  01  =  0. 

03  9  IB=IA 

040  DO  17  J*  1 , N 

041  IB*IB+1 

042  01*Ul4F  (J)  *CCKJS  (P  (IB) ) 

04  3  17  CONTINUE 

044  l^REA!  (U 1)  /  (LX*ETA*SIN  (?R)  ) 

045  WRIT 1(21, 18)  T 

046  18  FOBMAK*  TBANSNISSICN  COEFFICIENT  T-»,I14.7) 

C47  CT*DX*DY/(PI*ETA) 

C48  CT*CT*CT/(3.*PI) 

04  9  £TH*PI/(NTH-1) 

050  1BITE  (2  1 ,23) 


'•v-v-l 


v  v  v  *.*  •> 


*>  \»  *.«  *  «‘*  • v  *  *»  »'*  •  *  «  *.«.*»  .  V  V  >  '  . 


051  23 

FOBS AT (*  ANGLE *#4X#,TA01*#7I»,TA02,J 

052 

DO  19  J-1.NTH 

053 

lH-(J-l)  *DTH 

054 

CALL  PLANE  (IB.IX.DX ,DT. P) 

055 

TH*TB*P8 

056 

J1*0  • 

057 

CO  20  K*  1#  2 

058 

01*0. 

059 

DO  21  1-1.  N 

060 

J1-J1+1 

061 

OI-OUP(JI)  *V  (I) 

062  21 

CONTINOl 

063 

B*01*C0BJ<5  (01) 

064 

TAO  (X)  -CT*H 

06  5  20 

CONTINOl 

066 

MBITS  (21.22)  TH,  (TAO  (I)  .1*1,2) 

067  22 

FOBS  AT  (1X.F7.  2  ,2E1 1.  4) 

06  8  19 

CONTINOl 

C69 

STOP 

070 

END 

INPOT  DATA  IN  THE  FILE  NAOTZ3.DAT 


5  1  IS  0.5000000E-01  Q.5000000E-0 1  0. 27000C0E+03 


OUT  POT  Dm  III  THZ  FILE  NA0TZ4.DAT 


LI  LI  HTH  DZ  DI  TH 

5  1  19  0*5000 0008*01  0.50000008*01  0.27000008*03 

I 

*0.  15318*02*0.65268*01  0.66468*01-0. 64638-0  1  0. 13128*01*0.62748*01 
P 

*0.10008*01*0.14718*06*0.10008*01*0.29  418*06*0.10008*01*0.44128*06 
COBPFIC11MTS  V  Of  8AGHETIC  COB BE NT  EXPANSION  FUNCTIONS 
0.45118*02  0.59168*03  0.62388*02  0.81538*03  0.62388*02  0.81538*03 
0.45118*02  0.59168*03 
T8ANSHISSICN  COEFFICIENT  T*  0.  11412638*00 
ANGLE  T 101  TAU2 

0.00  0.00008*00  0.21868*02 
10.00  0.58858*04  0.21868*02 
20.00  0.23088*03  0.21888*02 
30.00  0.50168*03  0.21908*02 
40.00  0.84628*03  0.21938*02 
50.00  0.  12288*02  0.21968*02 
60.00  0.16028*02  0.21998*02 
70.00  0.  19188*02  0.22028*02 
80.00  0.21298*02  0.22038*02 
90.00  0.22048*02  0.22048*02 
100.00  0.21298*02  0.22038*02 
110.00  0.  19188*02  0.22028*02 
120.00  0.16028*02  0.21998*02 
130.00  0.12288*02  0.21968*02 
140.00  0.84628*03  0.21938*02 
150.00  0.50168*03  0.21908*02 
160.00  0.23088*03  0.21888*02 
170.00  0.58858*04  0.21868*02 
180.00  0.20888*15  0.21868*02 


IV.  THE  SUBROUTINE  YMAT 

The  subroutine  YMAT (LX,  DX,  DY,  Y)  stores  by  colunns  in  Y  the 
matrix  Y  that  appears  In  (55)  and  (56).  The  first  three  argu¬ 

ments  of  YMAT  are  Input  arguments.  The  aperture  is  LX  subsections 
long  In  the  x  direction  and  one  subsection  wide  in  the  y-directlon. 

DX  is  kAx  and  DY  is  kAy  where  k  is  the  wave  number.  Ax  is  the  subsec¬ 
tion  length  in  the  x  direction,  and  Ay  is  the  subsection  width  in 
the  y  direction. 

Minimum  allocations  are  given  by 

COMPLEX  TCOX+1),  TX(LX+1),  YXX(N) ,  Y (N*N) 
where  N  is  given  by  (54) . 

Inside  DO  loop  16,  Ic(I-l)  of  (23)  is  put  in  the  computer  program 
variable  TC(I+1) .  Here,  I  is  the  index  of  DO  loop  16.  Also  Inside  DO 
loop  16,  Ix(I-l)  of  (29)  is  put  in  the  computer  program  variable  TX(I+1). 
The  logic  Inside  DO  loop  16  is  best  understood  by  building  up  a  table  of 
variables  in  YMAT  versus  expressions  in  terms  of  variables  in  Section  II. 


Variables  in  YMAT 

Expressions  in  Section  II 

I 

i+1  where  i  appears  in  (23)  and  (29) 

YU 

Yu  of  (12) 

XU 

xu  of  (13) 

XL 

x£  of  (14) 

Rl 

r±  of  (22) 

U1 

Ux  of  (24) 

U2 

U2C2  of  (25)  and  (32) 

U3 

U3  of  (26) 

U4 

|  of  (27) 

TC(I+1) 


TX(I+1) 


2e  *  of  (23)  and  (29) 
r3  of  (39) 
of  (40) 

Axu  of  (41) 

Axl  °f  <42> 

k  of  (43) 
yu 

Cx  of  (31) 

C3  of  (33) 

C4  of  (34) 

I  (i)  of  (23)  with  1  -  1-1 

c 

\  Of  (35) 

X3  of  (37) 

X4  of  (38) 

I  (i)  of  (29)  with  1  -  1-1 


Taking  advantage  of  the  fact  that  I  (1)  Is  even  In  1,  line  46  stores 

c 

Ic(-1)  in  TC(1).  Taking  advantage  of  the  fact  that  1^(1)  Is  odd  in  i, 
line  47  stores  Ix(-1)  1°  TX(1). 

Inside  DO  loop  20,  lines  51  and  52  put  in  YXX(J-l)  the  square 
bracketed  term  on  the  right-hand  side  of  (9)  with  j-i  -  J-2.  Now,  we 


J-i  -0.1.2 . V2 


Because  Y^  is  even  in  (j-i),  (62)  becomes 


35^  Y4j  -  YXX(|j-i|+l),  ( j-i|  -  0,1,2 . Lx-2  ( 

Inside  nested  do  loops  23  and  21,  line  59  puts  Yjj  of  (63)  in  the 

computer  program  variable  Y(I  +  (J-1)*N)  where  N  is  given  by  (54). 


oo  re 
002 
003 
004 

cos 
006 
007 
008 
009 
010 
C11 
012 
013 
014 
015 
016 
017 
018 
019 
020 
021 
022 
C23 
02  4 
025 
026 
027 
028 
029 
030 
031 
032 
033 
034 
035 
036 
037 
033 
039 
04  0 
041 
042 
043 
044 

045  16 

C46 

047 

048 

049 

050 

051 

052 

053  20 


LISTING  OF  THE  SOBBCOTINE  THAT 
SOBBOOTINE  THAT  (1X,DX,DY,T) 

COflPLEX  U, 01, 02# 03, 04, EX,TC  (100)  ,TX(100) 
COBPLEX  YXX(  100)  ,1(1600) 

DX2*DX*DX 

B*IX-1 

0*(0.,1.) 

04*. 1666667*0 


10*. 5*DY 

054 

JY*0 

100*10*0 X 

055 

CO  23  J* 1, N 

Y02*YO*YO 

056 

DO  21  I*  1 , N 

YO 3*. 3333333*10 

05  7 

JY*JY* 1 

104*. 2 5* 0X2* YD 3* YO 

058 

K*IAES  (J— I) *1 

CO  16  I*  1  #  LX 

059 

Y(JY)*YXX(K) 

IP*I*1 

060  21 

CONTINOE 

XO* (I—. 5)  *DX 

061  23 

CONTIHOE 

X02*X0*X0 

06  2' 

BE10BN 

XL*X0-0X 

06  3 

END 

X12*U«XL 
Hi*  (1-1)  *DX 
H2*B 1*H1 
B0 1*1. -.5*8  2 

01=801*81*  (1.-. 1666667*82)  *0 
02*  (Hi— H01*O)  *XOD 
03*-.5-. 5*81*0 
IX*2.  *  (COS  (8 1)  -0*SI8  (8 1)  ) 

R7=IL2*Y02 

H8*X02*Y02 

H3*SQ8T|87) 

84*SQ8T(B8) 

XX 0*  10* ALOG  (  (10*84)  /XO) 

A XL* XL*1L0G (  (TO*B3) /ABS (XL) ) 

ATO*  YO*AIOG  ( (10*84)/  (XX*B3)  ) 

C1*AXO-iXL*AIO 

C3*Y03* ( XO*B4-XL*B3) *.1666667* (XD2*AX0-XL2*AXL*Y02*AY0) 
C4»I03* (XO*B8-XL*R7) 

1C  (IP)*  (C  1*0 1  *02 *C3*03*C4*0 4)  *EX 

IXO*XO*AXO 

AXL*X1*AXL 

XI*.  5*  (10*  (84-83)  ♦  AXO-  AXL) 

XJ*Y0*  {.  8333333B-1*  (B8*84-B7*83)  ♦.  125*  (X02*BU-XL2*33) ) 

1  ♦. 125* (X02*AX0-XL2*AXL) 

X4*B  1*Y0D*  (R2*  Y04) 

II  (IP)*  (X1*01*B1*02*X3*03*X4*04)  *EX/DX 

CONTINOE 

1C  (1)*TC  (3) 

IX  ( 1)  *-TX  (3) 

CO  20  J*  2  ,  L  X 
JM* J- 1 
J?*J*1 

YXX  ( Jfl)  »•  5*  (TC  (J)-TX(JP)  ♦  (J-.5)  *TC  (J?)  *TX(JH)- 
1  (J-3.5)  *1C  ( J«)  )  ♦  (TC  (JP)  -2.  *TC  (J)  *TC  ( JN)  )  /DX2 

CONTINUE 


V.  THE  SUBROUTINE  PLANE 


The  subroutine  PLANE(TH,  LX,  DX,  DY,  P)  stores  in  P(l)  to  P(N) 
the  quantities 


.  ,kAx  cos  6^~nc  2 

1  (P^nc)a„  -  sin  9inc  ( - - — jrr — )  eJkiAx  C0S  9 


,inc 


2AxAy  'li  '0y 


kAx  cos  6 


inc 


i*l ,2, . . . ,N 


(6 


Where  (64)  comes  from  (8)  and  N  is  (L  -1) .  The  aperture  is  L  sub- 

X  X 

sections  long  in  the  x  direction  and  one  subsection  wide  in  the  y 
direction.  In  (64),  k  is  the  wave  number.  Ax  is  the  subsection 
length  in  the  x  direction,  and  Ay  is  the  width  of  the  aperture  in  the 
y  direction.  Moreover,  the  subroutine  PLANE  stores  in  P(N+1)  to  P(2N) 
the  quantities 


1  /pin<\ 

2AxAy  1  i  xx 


.  .  ,inc 

.  kAy  cosj)  . 

_  8ln  2  >  j  (kAy/ 2)  cos  <|»lnc 

kAy  cos  4>lnc 


(6 


inc 

where  (65)  comes  from  (16).  The  angle  if)  in  (65)  is  assumed  to  be 

jLnc 

the  same  as  6  in  (64).  The  first  four  arguments  of  PLANE  are  input 

inc  inc 

arguments.  In  radians,  TH  is  the  common  value  of  9  and  $  in  (64) 
and  (65).  Furthermore,  LX  is  Lx>  DX  is  kAx,  and  DY  is  kAy.  LX  is  a 
positive  integer  greater  than  or  equal  to  2. 

The  minimum  allocation  for  P  is  given  by 

COMPLEX  P(2*N) 
where  N  is  given  by  (54) . 

Lines  5  and  6  set 
CX  -  kAx  cos  (TH) 

CY  ■  -J  kAy  cos  (TH) 
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If  cos  (TH)  ^  0,  then  lines  11  to  14  set 

sin  (“SJSSjffil) 

SX  -  sin  (TH)  ( - r-r - r - y 

kAx  cos  (TH) 


sin  (k^-c°S-lTH)) 

kAy  cos  (THY 
2 


If  cos  (TH)  *  0,  then  lines  8  and  9  set  SX  and  SY  equal  to  the  limits 
of  the  above  two  expressions  as  cos  (TH)  goes  to  zero,  namely 

SX  -  sin  (TH) 

SY  -  -  1 

Inside  DO  loop  13,  line  19  puts  in  P(I)  the  right-hand  side  of  (64)  with  i 
replaced  by  I  where  I  is  the  index  of  DO  loop  13.  Note  that  the  right- 
hand  side  of  (65)  does  not  depend  on  i.  Line  21  puts  in  U1  the  right-hand 
side  of  (65).  DO  loop  14  puts  the  right-hand  side  of  (65)  in  P(N+1)  to 
P(2*N)  where  N  is  given  by  (54). 


->V-' 


00  1C 
002 
003 
004 
005 
006 
007 
00810 
009 
010 
01 1  11 
012 
013 
014 
01  512 
016 
017 
018 
019 
02013 
021 
022 
023 
024 
02514 
026 


LZ STING  or  THE  SOB SOUTINE  PLANE 
SOBSCOTINE  PLANE  (TH,L!,DX,DY,P) 
COUPLE!  O,U1,P<1C0) 

CS-CCS  (TH) 

CX=EX*CS 
CX=.5*DY*CS 
IP  ICS)  11,10,11 
SX*SIH  (TH) 

ST3- 1. 

GO  TC  12 

SI*. S*CX 

SX— SIN  (SX) /SX 

SX=SIN  (TH)  *SX*SX 

SI*-SIH(CI)/CT 

U*{0.,  1.) 

N*LX-1 
CO  13  I-1.N 
S=I*CX 

P(I)  *SX*  (COS  (S)  ♦0*SIN(S)  ) 

01*SI*  (COS  (CT)  ♦  0*SIN  (CT) ) 

CO  14  J-1,N 
I»J  +  N 
P(I)  *01 
CONTINUE 
SETOHN 


THE  SUBROUTINES  DECOMP  AND  SOLVE 


The  subroutines  DEC0MP(N,  IPS,  UL)  and  SOLVE(N,  IPS,  UL,  B,  X) 
solve  a  system  of  N  linear  equations  in  N  unknowns.  The  input  to 
DECOMP  consists  of  N  and  the  N  by  N  matrix  of  coefficients  on  the 
left-hand  side  of  the  matrix  equation  stored  by  columns  in  UL.  The 
output  from  DECOMP  is  IPS  and  UL.  This  output  is  fed  into  SOLVE.  The 
rest  of  the  input  to  SOLVE  consists  of  N  and  the  column  of  coefficients 
on  the  right-hand  side  of  the  matrix  equation  stored  in  B.  SOLVE  puts 
the  solution  to  the  matrix  equation  in  X. 

Minimum  allocations  are  given  by 
COMPLEX  UL(N*N) 

DIMENSION  SCL(N),  IPS(N) 
in  DECOMP  and  by 

COMPLEX  UL(N*N),  B(N) ,  X(N) 


DIMENSION  IPS(N) 


in  SOLVE. 


More  detail  concerning  DECOMP  and  SOLVE  is  on  pages  46-49  of  [3]. 
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001  c 

LISTIHG  OF  THE  SOBBOOTINE  DECOBP 

002 

SOBBCOTIHE  DECCHP (H,IPS, OL) 

003 

COUPLE X  OL(1600) ,PIVCT,EH 

004 

EIHEHSIOH  SCL(40)  ,  IPS  (40) 

00  5 

DO  5  1*1,  H 

006 

IPS  (I)  »I 

007 

BH*0. 

008 

J1»I 

009 

DO  2  J*1,S 

010 

OLE3 AES  (BEAL  (UL(  Jl) ) )  +ABS  (AIBA6  (OL  (J1)  )) 

Oil 

J1*J1+H 

012 

IF (BH— OLH)  1,2,2 

013 

1 

BE3 OLE 

014 

A 

4 

COHTIHOE 

015 

SCL(I)*1./BH 

016 

c 

COHTIHOE 

017 

HHl*  H-1 

018 

K2*0 

019 

DO  17  R31,HH1 

020 

EIG30. 

021 

CO  11  I3K,H 

022 

IP3IPS  U) 

023 

IPK3IP*I2 

024 

SIZE3  (ABS  (BEAL  (OL  (IPK)  ))*KBS  (AIHAG  (OL  (IPK)  ) ) )  *SCL  (IP) 

025 

IF  (SIZE-BIG)  11,11,10 

026 

10 

BIG3SIZE 

027 

IP?3 1 

028 

11 

COHTIHOE 

029 

IF(IPV-K)  14,15,14 

030 

14 

J*IPS(K) 

031 

IPS  ( X)  3IPS(IPV) 

032 

IPS(IPV)3J 

033 

15 

KPP3IPS (K)  *K2 

034 

PIVOT3 OL  (KPP) 

03  5 

KP13K«-1 

036 

CO  16  I3K?1,H 

037 

KP3KPP 

038 

IP3IPS  (I) *K2 

03  9 

Efl3-0L  (IP) /PIVOT 

040 

18 

OL  (IP)  *-E8 

041 

DO  16  J»RP1,H 

042 

I?*IP+H 

C4  3 

RP*RP*H 

044 

OL  (IP)  »OL(IP)  ♦E3*0L  ( K?) 

045 

16 

COHTIHOE 

C4  6 

K2*K2+H 

04  7 

17 

COHTIHOE 

048 

BZTOBH 

045 

EHD 

.  •  -  <*>  vj*«  ,*•  ;.r*AVr-> 


050  C 

LIST IB G  OF  THE  SOBBCOTIBE  SCLVZ 

051 

SOBBCOTIBE  SOLVE  (S,IPS, OL, B#X) 

052 

CO  HP  LEX  OL  ( 1600)  ,B  (40)  ,X  (40)  ,S0H 

053 

DIBB1SIOB  IPS (40) 

054 

BP  1*1*  1 

055 

IP*IPS(1) 

056 

X  ( 1)  *B  (IP) 

057 

00  2  1*2,1 

058 

XP*IPS  (I) 

OS  9 

XPB*IP 

060 

IB  1*1-1 

061 

SO  (1*0. 

062 

CO  1  3*1, IB1 

063 

SDB*SOH+OL  (IP)  *X(J) 

064  i 

IP*IB+B 

C65  2 

1(1)  *B(IPB)-SOH 

066 

K2*B*fB-1) 

067 

IP*IPS  (B) *K2 

068 

X  (B)  *X  (H)  /OL  (IP) 

069 

00  4  IB ACS* 2, B 

070 

I*  BP  1— HACK 

071 

K2*K2— H 

072 

IPI*IPS (I)  *K2 

073 

IP  1*1*1 

074 

SOH*0. 

075 

IP-IEI 

076 

CO  3  J*IP1,» 

077 

IP*IP*B 

078  i 

SOB*SOH«OL(IF)  *X  (J) 

079  4 

X  (I)  *  (X  (I)  -SOB)  /OL  (IPI) 

08  0 

8  BTOBB 

081 

END 
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